Marked influence of the nature of chemical bond on CP-violating signature in 

molecular ions HBr^ and HI"*^ 
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Heavy polar molecules offer a great sensitivity to the electron Electric Dipole Moment (EDM). 
To guide emerging searches for EDMs with molecular ions, we estimate the EDM-induced energy 
corrections for hydrogen halide ions HBr+ and HI+ in their respective ground X^^/i states. We 
find that the energy corrections due to EDM for the two ions differ by an unexpectedly large factor 
of fifteen. We demonstrate that a major part of this enhancement is due to a dissimilarity in the 
nature of the chemical bond for the two ions: the bond that is nearly of ionic character in HBr"*" 
exhibits predominantly covalent nature in We conclude that because of this enhancement the 

HI'*' ion may be a potentially competitive candidate for the EDM search. 
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A non-vanishing permanent electric dipole moment 
(EDM) of a particle simultaneously violates two discrete 
symmetries: parity (P) and time reversal (T). By the 
virtue of the CPT theorem, the T-violation would im- 
ply CP-violation 0. Searches for EDMs of atoms and 
molecules provide important constraints on competing 
extensions to the standard model of elementary particles. 
For example, the most stringent limits on electron EDM 
come from a table-top experiment with atomic Tl 
As with atoms, the internal energy states of heavy polar 
molecules can show evidence of EDMs of the constituents. 
Compared to atomic experiments, where application of 
strong external E-field is required to enhance sensitivity, 
the experiments with polar molecules rather rely on the 
inner electric molecular field Si^t exerted upon the heav- 
ier atom. This field can be several orders of magnitude 
larger than the attainable laboratory fields. This notion, 
first elucidated by Sandars Q, has been exploited in ex- 
periments with YbF Q and TIF (E ^ ng] molecules. We 
also mention ongoing experiment with metastablc PbO 
molecule 0- 

A relatively small laboratory field is still required in 
the EDM experiments to polarize the molecule. Since 
the E-field would accelerate a charged particle out of an 
apparatus, EDM experiments are typically carried out 
using neutrals. It has been recently realized by Stutz 
and Cornell ^3 that this limitation may be overcome 
with ion traps: electrostatic force exerted upon the ion 
by the polarizing E-field can average to zero if the po- 
larizing field rotates rapidly in space, with the requisite 
spectroscopy then being performed in a rotating frame of 
reference. Moreover, the long coherence times in the trap 
would improve statistics compared to traditional beam 
and gas cell approaches. Because of this improved statis- 
tics, molecular ions with a relatively weak sensitivity to 
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electron EDM could provide competitive constrains. In 
particular, the hydrogen halide ions HBr"^ and HI+ m 
their lowest rovibrational state of the ground 
term are considered as attractive candidates for the pro- 
posed experiment pdil |. 

The goal of this paper is two-fold. Firstly, we provide 
a guidance to emerging EDM searches with molecular 
ions ^3 by computing EDM-induced energy corrections. 
Secondly, we elucidate the important role of the chemi- 
cal bond in enhancement of electron EDM in molecular 
systems. While both HBr'*' and HI^ ions have a simi- 
lar electronic structure, the chemical bond in HBr'*' is of 
ionic nature, while for heavier HI'*' it is predominantly 
covalent [ill . We find that this evolution in the char- 
acter of the chemical bond has a marked effect on the 
EDM-induced energy corrections. From the experimen- 
tal point of view, our computed EDM-induced energy 
correction for HBr+ is too small to produce competitive 
bounds on the electron EDM in experiment [lOi] . By con- 
trast, the pronounced covalent bond enhancement for the 
111+ ion, illuminated here, makes it a potentially compet- 
itive candidate for the electron EDM search. 

Molecular structure and EDM-induced corrections. 
The molecular structure of low-lying rotational states of 
hydrogen halide ions HBr'*' and HI+ can be well clas- 
sified by the Hund's case (a). Relativistic effects split 
the ground X ^11 electronic term into two components: 
^113/2 and distinguished by fi, projection of the 

total electronic angular momentum along the molecular 
axis. ^113/2 is the ground electronic term and it is con- 
sidered as a possible candidate for the EDM experiment. 
In the estimates below we will employ the following val- 
ues of the equilibrium internuclear separations [ ill Il2| : 
Re « 1.448 A for HBr+ and Re « 1.632 A for HI"^. Un- 
less noted otherwise, atomic units, h = |e| — me = 1, are 
used throughout the paper. 

In Hund's case (a) the molecular eigenfunctions 
including the nuclear rotation can be described as 
\AT,n; JMj) = \ JMjVl) IAS 17) , where A and S are pro- 
jections of the electronic orbital momentum and spin onto 
the internuclear axis, J is the total molecular momentum 
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(including nuclear rotation) and Mj is the laboratory 
frame projection of J. The rotational part | JMjfi) may 
be expressed in terms of the Wigner D functions. While 
in the lowest-order JMj) and |-A - S - fi; JMj) 

states have the same energies, at the finer level each ro- 
tational state of the ^Iln terms splits into so-called A~ 
doublet |l3j| due to rotational and spin-orbit perturba- 
tions. The eigenstates of the field-free molecular Hamil- 
tonian (disregarding EDM) are e// parity states, com- 
posed of linear combinations of the two above states. 

An externally applied electric field £q couples the e// 
parity states. For a sufficiently strong E-field the eigen- 
states can be classified by a definite value of 57, rather 
than by the e// parity label. In this case the correction 
to the energy due to electron EDM can be parameterized 
as Q 5W{ J, M, n) = Wd^- It is defined as an expecta- 
tion value 

Wdfl = {AT.n;JMj\He jASll; JMj) = (A S 17|iJe|A E rj) 

(1) 

Here iJg = —de{jo — 1)S • £ is the pseudo-scalar cou- 
pling 13 of an electron EDM de to an electric field £ 
(this internal molecular field is to be distinguished from 
the externally applied field). The expectation value ^ is 
accumulated in the region of strong fields, i.e., mainly in 
the vicinity of the nucleus of the heavier halogen atom. 
A common approximation is that the electric field is 
produced by a spherically symmetric charge distribution 
£{r) « Z/r^r, where Z is the nuclear charge of the 
heavier atom, and r = coincides with its center. 

Chemical bond. In the following we make an order-of- 
magnitude estimate of the EDM factor Wd using a qual- 
itative model of an isolated atomic particle perturbed by 
its molecular counterpart. In this regard it is important 
to discuss the nature of the chemical bond in the hy- 
drogen halide IIX+ ions. It can be be described by two 
limiting cases 12]: ionic (H"'": X) and covalent (H • • ■ X"*") 
bonds. In the case of the ionic bond the halogen atom is 
electrostatically perturbed by a proton. When the bond 
is covalent, the halogen atom is singly ionized (■^P state), 
while the hydrogen atom is in its ground state. 

Although both HBr+ and III+ ions dissociate to the 
covalent limit, the chemical bond at intermediate separa- 
tions can be better characterized from molecular spectra. 
In particular, the hyperfine structure is of significance to 
our consideration, because both the EDM coupling and 
the hyperfine interaction are sensitive to behavior of the 
molecular orbitals near the nuclei. An analysis of the hy- 
perfine structure in Ref. indicates that the bond for 
the IIBr+ ion can be adequately described as being of the 
ionic nature. As to the III+ i on, the hyperfine-structure 
analysis by the same authors [ij shows that the bond is 
predominantly of the covalent character. 

Below we consider both ionic and covalent bonds. Our 
semi-qualitative calculations follow a general scheme sim- 
ilar to those described in Ref. 0, 0|- Firstly, we de- 
termine the effective molecular electric field fmt exerted 
upon the heavier halogen atom/ion. Then we use the 
first-order perturbation theory in the interaction with 



fint to determine mixing of the atomic states of oppo- 
site parity. Finally we compute the expectation value 
of the EDM-coupling operator using ab initio relativistic 
atomic structure codes. 

Ionic bond approximation for the HBr^ ion. In the 
case of the ionic bond the halogen atom is electro- 
statically perturbed by a proton. In the spirit of the 
LCAO 19] method we expand the electronic wavefunc- 
tion in terms of atomic states $i of the halogen atom 



(2) 



where the total angular momentum J^^i of the atomic 
state $i and its projection on the molecular axis Me_i 
are constrained to Je,i > \^\ and M^^i = ri. We deter- 
mine the expansion coefficients Ci using the first order 
perturbation theory in the interaction V due to electro- 
static field exerted upon the halogen atom by the proton 
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where $0 is the ground atomic state of the proper sym- 
metry and Ei are the energies of atomic states. 

Keeping only the leading dipole term in the multipole 
expansion of the interaction of atomic electrons of the 
halogen atom with the proton, the perturbation V w 
— D • £iut, £int being the electric field of the proton at 
the position of the atom and D the atomic electric dipole 
operator. It is this strong electric field that produces 
enhancement of the electron EDM in molecular ions. 

Finally, the EDM-induced energy correction is 
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In the following we will use a shorthand notation 



T ^ He {Eq — Ha) 



with Ha being the atomic Hamiltonian so that 
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It is worth noting that all the quantities (except for 
empirical Re) in the Eq. Q are atomic ones and we 
employ atomic-structure methods to evaluate this sum. 
First we employ Dirac-Hartree-Fock (DHF) approxima- 
tion and then more elaborate configuration-interaction 
(CI) method. All calculations carried out here arc ab 
initio relativistic. 

The halogen atoms Br and I are open-shell systems 
with one hole in the outer np3/2-shell, n = 4 for bromine 
and n — 5 for iodine. In the DHF approximation 
the atomic orbitals |i} satisfy the eigenvalue equation 



3 



^dhfK) = Eil*), where the Dirac Hamiltonian /idhf in- 
cludes an interaction with the field of the nucleus and the 
self-consistent field of the electrons. In the DHF approx- 
imation we obtain 

^^^^2^{jMm^^ (7) 

e j ^9 

where g denotes the hole state and the summation 

over i extends over a complete set of orbitals, including 
both core and virtual orbitals. 

Numerically we carried out the summation using the 
B-spline pseudo-spectrum technique 16J. The pseudo- 
spectrum was generated using the DHF potential of the 
ground ^^3/2 atomic state. In a typical calculation we 
used a set of basis functions expanded over 100 B-splines, 
which provided numerical accuracy sufficient for the goals 
of this paper. Among other technical details it is worth 
mentioning that while integrating the radial Dirac equa- 
tion, we used the potential produced by a nucleus of the 
finite size. 

To investigate a potentially large correlation effects be- 
yond the DHF approximation, we have also carried out 
configuration-interaction (CI) calculations for Br within 
the active space of seven 4s^ Ap^ valence electrons. In this 
method, the many-electron wave functions were obtained 
as linear combinations of determinants composed from 
as single and double excitations of the valence electrons 
from the active space. Finally, following Dalgarno-Lewis- 
Stcinhcimer method |l7j . we carried out the summation 
over intermediate states in Eq. (|12|l by solving inhomoge- 
neous many-body Dirac equation inherent to the method 
and computed the sum Q . More details will be provided 
elsewhere. 

The resulting DHF value of the EDM-induced energy 
correction (hole in the 4p3/2 shell ) 

t^df^[HBr+, ionic, DHF, X ^3/2] = -1.5 x lO^^d^ , (8) 

A similar DHF calculation assuming a hole in the 'Api/2 
shell leads to a lOO-fold increase in the value of the EDM 
correction 

H/dl7[HBr+, ionic, DHF, X ^1/2] = 1.64 • (9) 

A large difference in the values of the Wjil parame- 
ter for the two cases can be explained as follows. The 
EDM-couphng operator Hf. is a pseudo-scalar: it does 
not change the total angular momentum of a state, but 
fiips its parity. For example, if the hole state g has Pzj2 
angular character, then the intermediate states in Eq.jTl) 
are ^3/2 orbitals. Similarly, the Pi/2 hole state requires 
Si/2 intermediate states. It is well known 2j, that since 
the states of lower orbital momentum have a larger prob- 
ability to be found close to the nucleus, this selection rule 
has a profound effect on the order of magnitude of the 
EDM factor Wa- 

One may argue that an enhancement of the EDM fac- 
tor for the X^l3/2 state may arise due to particle- hole 



excitations, when si/2 (pi/2) electron is excited from the 
core to the pi/2 (si/2) orbital. It is easy to demonstrate 
in the DHF approximation, that while the individual con- 
tributions from such excitations are certainly large, their 
sum vanishes. It is the reason why the closed-shell sys- 
tems are largely insensitive to the electron EDM Q, i.e 
the EDM-induced energy correction arises only due to an 
unpaired electron. Correlations (many-body effects be- 
yond DHF) may potentially spoil the presented argument 
and we have carried out the correlated CI calculations. 
The result, 

I^df^ [HBr+, ionic, CI, X ^3/2] = -2.6 x lO^^de , (10) 

is of the same order as the DHF value. 

As a reference, here we also present the DHF value for 
the HI+ ion in the ionic bond approximation 

H/dl7[HI+, ionic, DHF, X 2H3/2] = -7.0 x lO^^^ . (n) 

Covalent bond approximation for HI+ . From the pre- 
ceding discussion it is clear that a participation of the 
unpaired Pi/2 or Si/2 orbital in the ground-state configu- 
ration of the heavier molecular constituent is important 
for gaining large Wd parameter. Qualitatively we can 
hope that such an enhancement for the X^^^2 compo- 
nent may arise when the chemical bond acquires cova- 
lent character (case of HI+ ion). Indeed, in the covalent 
bond approximation, the halogen atom becomes singly 
ionized its ground state being ^P. The ground state has 
two p holes in the outer shell, so that the corresponding 
relativistic many-body states are composed from linear 
combination of p^y^2-Pi/2' P3/2P1/2 ^^"-^ P3/2P3/2 single- 
electron configurations (the superscript —1 designates a 
hole state). Therefore the unpaired j — 1/2 orbital be- 
comes involved in the calculations, and indeed, as shown 
below, this leads to a significantly larger EDM-induced 
energy correction for the X^^^2 term. 

The HI"*" ion may be pictured as the iodine ion 1+ in the 
•^P state perturbed by the neutral hydrogen atom in its 
ground state. First let us derive the internal electric field 
fint and the associated mixing of opposite parity states 
of the iodine ion. Qualitatively, the field of 1+ induces a 
dipole moment of the hydrogen atom \Dh\ = ao/Rl, 
where ao = 9/2 is the polarizability of the hydrogen 
ground state. In turn, the induced dipole moment exerts 
a field at the position of the iodine ion £int — 2ao/R^z. 
Thus the iodine ion is perturbed hy V ^ —2aoDz/Rl , 
where D is the atomic dipole moment operator for I"*". 

Again we limit our consideration to a qualitative esti- 
mate and use the first-order perturbation theory in the 
molecular field, so that the EDM-induced energy correc- 
tion is 

W^n = ($o|T|<i>o) , (12) 

where the operator T is given by Eq. , except now all 
the participating operators in that expression are to be 
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understood as being for the iodine ion and $o is its prop- 
erly symmetrized ground state. The above expression 
differs from the analogous formula ^ for the ionic bond 
by a prefactor characterizing the internal molecular field 
Sint acting upon the halogen atom/ion. Compared to the 
ionic bond, this perturbing field becomes 70% weaker. 

In order to carry out the calculations with the rela- 
tivistic operator He, we express the unperturbed non- 
relativistic molecular wavefunction in terms of the rela- 
tivistic wave-functions of the iodine ion, |^Pj, and 
the hydrogen, |1si/2,Mh), 
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Since the expectation value of the EDM coupling oper- 
ator He is accumulated close to the nucleus of the heavy 
iodine ion, H^ is essentially a one-center operator and a 
generalization of Eq. (|12|l for the two-center wavefunction 
((Till reads 
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m 



"113/2 |^^ePn3/2)covalont 
1 



(14) 



- (3P2, 1| T |3F2, 1) + - f Pi, 1| T |3Pi, 1) + 
i (3P2, 1| T 1) + 1 (3Pi, 1| T |3P2, 1) , 

We calculated the values of matrix elements for 
iodine ion within the CI approach similar to the 
one described above for Br. The computed values 
are (^Ps, l| T |3p2, l) = 6.44, (^Pi, l| P pPi, l) = 
-13.44, (3Pi,l|T|3p2,l) = -11.24, and 
(3P2,l|PpPi,l) = 2.54- Finally, 



H^dJ^[HI+,covalent,CI,X2n3/2] = -0.4d, 



(15) 



We notice a sizable enhancement compared to the value 
of — 7x 10~^4 obtained in the ionic bond approximation. 

Conclusions. First of all the EDM-induced energy cor- 
rection for the X^l3/2 state of HI+ is about 15 times 
larger than for IIBr+. A lesser part of this enhancement 
comes from the well-known scaling of CP-violating 
matrix elements 2], when bromine [Z — 35) is replaced 
by the heavier iodine {Z = 53). A more substantial fac- 
tor, illuminated in this work, is the evolution in the na- 
ture of the chemical bond. To reiterate, the CP-violating 
matrix elements are much larger for the pi/2 states than 
for P3/2, due to the fact that the values of the relevant 
matrix elements are accumulated close to the nucleus. In 
the ionic bond case of HBr"*" , the EDM correction arises 
from an unpaired P3/2 hole state in the outer shell and 
CP-violating effects are suppressed. By contrast, the co- 
valent bond of III+ in addition opens the pi^2 shell, lead- 
ing to a marked enhancement. 

Typical values [THIl of the EDM-induced energy cor- 
rections for heavy neutral polar molecules PbO and YbF 
are on the order of 104 atomic units. Our computed 
value for III+ is an order of magnitude smaller. Yet, 
when compared with the conventional beam and gas -cell 
experiments, the proposed trapping experiment \UJ^ has 
a better statistical sensitivity so that molecular ions with 
smaller enhancement parameters, such as HI"*", may suf- 
fice. By contrast, the EDM correction for HBr+ is too 
small to be of experimental interest. As shown here, it 
is the covalent bond of HI+ that makes this ion a poten- 
tially competitive candidate for the emerging searches for 
EDMs with molecular ions. 
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